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O . Abstract 

o 

CN | We use the Omnes representation to obtain the q 2 dependence of the form factors 

■ f +,0 (q 2 ) for semileptonic H — > it decays from elastic irH — > irH scattering amplitudes, 

where H denotes a B or D meson. The amplitudes used satisfy elastic unitarity and are 
obtained from two-particle irreducible amplitudes calculated in tree-level heavy meson 
chiral perturbation theory (HMChPT). The (^-dependences for the form factors agree 
with lattice QCD results when the HMChPT coupling constant, g, takes values smaller 
than 0.32, and confirm the milder dependence of /° on q 2 found in sumrule calculations. 

m 

<N ; 1 Introduction 

In this letter we present a description of the form factors / + and / describing semilep- 
tonic H — > Ti decays, where H denotes a D or B meson. For the B meson this exclusive 
semileptonic decay can be used to determine the magnitude of the CKM matrix element V u b, 
currently the least well-known entry in the CKM matrix. Ultimately, experimental mea- 
surements of fg{q 2 ) for given momentum-transfer q will be compared directly to theoretical 
determinations at the same q 2 values to determine \V u b\. In the interim, it may be helpful to 
consider the decay rate integrated partially or completely over q 2 , but this requires knowl- 
edge of the q 2 dependence of the form factors. Lattice calculations and sumrule calculations 
apply in (different) restricted ranges of q 2 while dispersion relations may be used to bound 
the form factors over the whole q 2 range |L], [|, or as a basis for models ||. A variety of 
models exists for the whole range of q 2 . One can ensure that general kinematic relations 
and the demands of heavy quark symmetry (HQS) are satisfied, but an ansatz, such as pole, 
dipole or other forms, is still required [|], [5]. 

Here we use the Omnes representation to obtain the full q 2 dependence of these form 
factors from the elastic ttH — > ttH scattering amplitudes. For our application we have an 
isospin-1/2 channel, with angular momentum J = 1 or for / + and f° respectively. We 
rely on the following description of the (inverse) amplitude for elastic nH — > nH scattering 
in the isospin /, angular momentum J, channel, with centre-of-mass squared-energy s and 
masses m and M respectively 0, 

TfJis) = -7 (a) - Cjj + l/V u (s), (1) 

where Vu is the two-particle irreducible scattering amplitude and Cjj is a constant. Cu 
and Vu are real in the scattering region. This description implements elastic unitarity auto- 
matically. Equation (|l|) is justified by a dispersion relation for T _1 , where the contributions 
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of the left hand cut and the poles (if any) are contained in — Cjj + 1/V/j. io gives the exact 
contribution from the right hand cut, after any necessary subtractions^ The description of 
equation (JJ) may also be justified by an approach using the Bethe-Salpeter equation. 

Once Tjj is known, we can compute the corresponding phase shift 8ij. In turn, 5u can be 
used in an Omnes representation || giving fu(q 2 )/ fu(0) in terms of an integral involving 
the phase shift, assuming that at threshold the phase shift should be nn, where n is the 
number of bound states in the particular channel considered, and Su(oo) = kn, where k 
is the number of zeros of the scattering amplitude on the physical sheet (this is Levinson's 
theorem 0). 

We determine Vu from tree level heavy meson chiral perturbation theory (HMChPT) 



TP] , which implements HQS and is a double expansion in powers of 1/M and momenta, 
where M is the heavy meson mass. The parameter Cjj in equation (JJ) partially accounts 
for higher order contributions in the expansion ||. 



We find consistency of our description with lattice results for the D — > n [fLl|, [12] and 



B — > 7r WA O] form factors if we set the HMChPT coupling, g, to values smaller than 0.32. 



This upper bound is in reasonable agreement with other determinations, but g is not very 
well known IH3L 



Our model and the Omnes representation are not guaranteed at high energies where 
inelasticities become important. However, our hypothesis is that only the low-lying states 
and energies should influence the form factors we consider. 

A dispersive approach to the / + form factor was taken by Burdman and Kambor H, who 
also used HMChPT to calculate the phase shift in nH — ► nH scattering. Here by working 
with the inverse amplitude we can ensure that Watson's theorem and elastic unitarity are 
satisfied exactly. Moreover, we compute / + and /° together to examine whether different be- 
haviours in q 2 are found, consistent with lattice QCD results and allowing extra information 
from f° to be used to constrain / + . 



2 Scattering Amplitudes and Form Factors 

We compute Vi/ 2 , the two-particle irreducible amplitude for irH scattering in the isospin 
1/2 channel, 7r(pi)H(Mv) — ► Tr(p 2 )H(Mv + q 2 ). Here, v is the four-velocity of the initial 
heavy meson of mass M. The pion mass is m. We use the direct tree level interaction from 
the lowest order HMChPT lagrangian, together with tree diagrams for H* exchange which 
involve the leading interaction term with coupling g [TO, ITJ| . The result is, 



V1/2 = +V-p 2 + g 2 (pvP2 - V-pxV-p 2 ) ( 3 _ . + ~T~a" J f " ( 2 ) 



Here, / = 130.7 MeV is the pion decay constant and A = (M 2 - M 2 )/2M % .\/, - 
where is the heavy vector meson mass. We subsequently project V1/2 onto the angular 
momentum and 1 channels. 

The full scattering amplitude at centre of mass energy-squared s, in the isospin I and 
angular momentum J channel, is obtained in our approach from equation (ffl). The phase 



1 Iq is calculated from a one-loop 'bubble' diagram. In the notation of reference M, Io(s) = Ta((r 
M) 2 ) — Tg(s), where M and m are the masses of the two propagating particles 
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shift 5jj is then obtained from, 



where X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx) is the usual kinematic function. 

Once the phase shift is known, we use the Omnes representation to obtain the q 2 depen- 
dence of the form factors as follows: 



H'j 2 ) 

7m =exp 



q 2 r°° 5u(s)ds 



IT J(rn+M) 2 s(-S — q 2 ) 



(4) 



In this work, we always have / = 1/2. The form factor / + is obtained when J = 1 and 
depends on J p = 1~ resonances, while /° is obtained when J = and depends on J p = + 
resonances. We perform the integral numerically, taking the upper limit as 100 times the 
lower limitf]. The form factors are equal at q 2 = 0: / + (0) = /°(0). 

• For J p = 1~ we take C = for the D decay because the D* resonance is so close 
to threshold that we expect it to saturate all the counterterms in HMChPT (compare 
to vector meson dominance in irn scattering in ordinary chiral perturbation theory). 
Calculating C in this case reveals the value C = 8 x 10~ 6 . We still have the freedom 
to vary the lowest order coupling constant g in HMChPT. For the B meson decay, we 
set C = —0.0014 to keep the B* pole at its correct mass. 

• For J p = + we ignore D* and B* s-channel exchanges, which have the wrong quantum 
numbers to contribute in this case. These exchanges only contribute because of the 
heavy meson mass expansion implicit in HMChPT. Instead we keep C non-zero, setting 
C = —0.0051 for the D-physics case to get a resonance at about 2350 MeV, and C = 



-0.0016 for S-physics to get a resonance at about 5660 MeV [16|. 



The values of C are determined by demanding that T" 1 (ReT -1 ) vanishes at the position 
of a pole (resonance). For the J = 1 channels, V^ 1 vanishes by construction at the positions 
of the D* or B*, and so, from equation ([[]), C is independent of g. In the J = channels, 
(7-dependence enters in V" 1 , but only through the t-channel tree graphs, and is very weak. 
C varies by less than 0.5% for < g < 0.45 in the D-meson case and the dependence is even 
weaker for the 5-meson case. 

We noted that in using the Omnes representation || of equation @, the phase shift 
at threshold should be mr, where n is the number of bound states in the channel under 
consideration. Thus n = in all channels used here except for J p = 1~ in the B case where 
n = 1 to account for the B*. In fact, our model also gives a bound state in the + channel in 
the B case, which we ignore. One could try to improve the model to avoid this unphysical 
bound state by replacing C with a function of q 2 (the function should have no right hand 
cut). 



2 /r>7r(9max)j where q^ax — i m D — rn-Tr) 2 , varies by less than 1% as the upper limit of integration varies 
from 50 to 200 times the lower limit, and the variation is smaller at lower q 2 . 
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Figure 1 Phase 5dk f° r the I J = 1/2,1 (left) and 1/2,0 (right) channels in Dtt scattering. The 
inset on the left shows the resonance at y/s = mp* = 2010 MeV. Phases are calculated with 
5 = 0.21. 



3 Semileptonic Decays 

The process D* — > Dtt is kinematically allowed, so the D* is a resonance in Dtt scattering. 
In HMChPT the decay rates of D* + to D°it + and D + ir° are given to lowest order by, 

T{D*+^D\+) = ^ T(D*+ -> D+vr ) = (5) 

071"/ z lZ7TJ z 

The sum of these rates can also be obtained from the slope of the phase shift at the resonance 
mass. We find that these two methods agree for a range of g values. 

The D* exchange is included in our tree level amplitude, and we expect it to saturate the 
counterterms in HMChPT, so in calculating Ti/2,1 we set C = as noted above. Figure [l] 
(left) shows the phase shift obtained for J = 1. With input masses, mp = 1864.5 MeV and 
tud* = 2010 MeV, the D* resonance shows up as the jump of it in the phase at the D* mass. 
In the J = channel, we tune C to produce a resonance at the expected mass of the Dq at 
2350 MeV |I6| . The J = phase shift is shown on the right in figure |1|. 



In the B case, the decay process B* — > Bit is not kinematically allowed and the B* 
meson is a pole, sitting between the maximum physical q 2 value for the form factor, g^ ax = 
(wb — m^) 2 , and the start of the physical cut at q 2 = (mB + m-n) 2 . Again we use the physical 
pseudoscalar and vector meson masses as inputs, vtlb = 5278.9 MeV, uib* = 5324.8 MeV. 
The phase shift for the J = 1 case is shown on the left in figure 0. The appearance of the B* 
as a bound state between and q 2 = {m B + m n ) 2 is signalled by the vanishing of 1 (s) 
at y/s = m B *. The J = phase shift appears on the right of figure 0. 

From the phase shifts we find the form factors / + and /°. We perform a simultaneous 
three-parameter fit to the UKQCD and APE lattice results JTT], [12], [13| for the form factors 
f + (q 2 ) and f°{q 2 ) which determine the B and D semileptonic decays. The free parameters 
are the HMChPT coupling constant g and the form factors at q 2 = 0: /s(0) for B — ► 71 
decays and /d(0) for D — » 71 decaysQ. The best fit parameters with 39 degrees of freedom 
are 

g = 0.21 + §;^}, / B (0) = 0.39 ± 0.02, f D (0) = 0.60 ± 0.02 with x7 d o f = 0.34. (6) 



3 To use the results in (II), we take Zy^ = 0.88 for the vector renormalisation constant connecting lattice 
and continuum results. 
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Figure 2 Phase 5b-k for the IJ = 1/2, 1 (left) and I, J = 1/2,0 (right) channels in Btt scattering. 
Phases are calculated with g = 0.21. 
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Figure 3 Form factors in D —* ir (left) and B — > tt (right) semileptonic decays. The squares 
(circles) denote / + (/ ) from lattice calculations, while the long-dashed (short-dashed) lines denote 
the fitted curves for / + (/ ). Solid symbols are results from UKQCD O, [l3|], open symbols are 
results from APE 11T 



Results can be seen in figure |^. Errors in the fitted parameters are statistical and have been 
obtained by increasing the value of the total \ 2 by one unit. A word of caution must be 
stated about the results for the HMChPT coupling constant g. Scalar channels are almost 
insensitive to this parameter. For the vector channels, in the case of D meson decay, the 
resonance is so close to threshold that it completely dominates the process, independent of 
the value of g, as long as the resonant contribution is more important than the background. 
This turns out to be true as long as g is greater than 0.001, thus the smallest value g can 
take is 0.001 and not zero as can be inferred from equation (|6]). To clarify the dependence 
of our results on g, we show in figure f| both \ 2 an d /b(0), /d(0) versus g, for g > 0.001. 
In the first figure the line at \ 2 = 13.28 shows the minimum value of x 2 ; while the line at 
X 2 = 14.28 determines the upper error. We also show best fit values, with fixed g, of /s(0) 
and /d(0) versus g. The points with errors correspond to the results quoted in equation (|6|). 

We note that f£ is well-approximated by a simple pole form with the D* giving the pole 
mass, while /£, is noticeably 'flatter' in q 2 . This is consistent with lattice results. For the 
B — >• 7i case, /g is well-approximated by a pole form with the pole mass of order the B* 
meson mass. The f% form factor has much less q 2 dependence, consistent with the behaviour 
found in lattice calculations. 



5 



X 




/(O) 




0.1 0.2 0.3 0.4 0.5 0.6 

9 



Figure 4 Left: chi-squared for the fit described in the text as a function of the HMChPT coupling 
g, for g > 0.001. Right: values of /g'°(0) (lower curve) and /d'°(0) (upper curve) as functions of 
g. The points with errors on the right are the best fit values of equation (0) at g = 0.21. 



We have also determined the coupling g and form factors at q 2 = separately for D and 
B decays using independent fits to the UKQCD and APE lattice data for D and B. The best 
fit values turn out to be the same as in equation (H), although gr> can be as large as 0.46 while 
still giving an acceptable chi-squared. To compare with light cone sumrule (LCSR) results, 
we take the LCSR values f d*Qd* d-k = 2.7 ± 0.8 GeV and f b*Qb*b-k = 4.4 ± 1.3 GeV [|i~7| . 
and combine with lattice calculations of the vector meson decay constants from Becirevic et 
al and UKQCD |L9|, to yield 



9d 



o.35 + 8;l? 

0.39 ±0.12 



f D * from [18] 
f D . from [0| 



9b 



0.23 ±0.08 

n oc ± 0.10 
- 0.09 



f B * from |T§ 
fs* from |19| 



(7) 



The values are quite compatible in the B case, less so for D decays, although, as noted above, 
our fit for allowed a large variation above the best-fit value. The value of /d(0) found 
here agrees well with the LCSR result /-(0) = 0.65 ±0.11 |I7| , while /b(0) in equation (|) 
is higher than the LCSR value /- (0) = 0.28 ± 0.05 [ 17| . In the D case, the D* resonance is 
only a few MeV above threshold and the range of q 2 for the semileptonic decay is not large, 
so one expects a simple pole form for / + to work well. For B physics, the effects of higher 
resonances and continuum states are evidently more important: such effects are incorporated 
in LCSR calculations but are not present in the very simple model used here. We address 
this issue in section |] below. 

Heavy quark symmetry (HQS) is an input in HMChPT. The HQS scaling relations for 
the B decay form factors at g max are preserved because / + (9 max )//°(9max) i s proportional to 



exp 



y 

imax 
7T 



5 + -5 c 

2 

max/ 



(M+m)2 S( S - g max ) 



ds 



The above result relies on the equality of the form factors at q 2 = 0, / + (0) = /°(0). If 
5 + — 5° = 7T, which we see is satisfied by our phase shifts at large i/i, then the ratio 
/ + (o , max)//°(5 , max) = V(l — 9max/(^ + m ) 2 ) as demanded by HQS, where M and m are the 
masses of the heavy meson and the pion respectively. 

We have applied the same approach to describe semileptonic D — > K decays. Here, 
it gives form factors flatter than lattice results |Ll]| and the experimental evidence p0 . 
However, corrections of both types mx/mo and m 2 K / (Air f n ) 2 to the tree level HMChPT 
results used here are expected to be sizeable in this case. 
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Figure 5 Variation of fg (0) with g and. (/res 

where the B* coupling is g and a second resonance is 
added in the J = 1 channel with coupling g ies . On the left the resonance mass is m res = 6100 MeV 
and contours are plotted from /g(0) = 0.20 to 0.28 in increments of 0.01; on the right m res = 
8100 MeV and contours are plotted from /^(0) = 0.28 to 0.37 in increments of 0.01. 



4 Extra Resonances 

We noted above that our result for /b(0) in equation (||) is higher than the LCSR value 
of around 0.28, while our fit for fsiQ 2 ) is well-approximated by a simple pole form with 
pole mass of order mg.. This suggests that deviations from B* pole dominance can become 
significant at low q 2 . This phenomenon was also noted by Burdman and Kambor H who 
implemented a constrained dispersive model for f^. Likewise, lattice results have favoured 
dipole forms in fits to f£ § [L| [T3|. 

To address this issue we have added a second resonance of mass m res by hand in the 2PI 
J = 1 amplitude Vi/2,1, coupling it like the H* but with its own coupling strength g res . In 
the D-meson case, we already had a good fit to the lattice results and a consistent value for 
/d(0). If m res is large enough the extra resonance does not disturb this picture. In the B 
case, we can easily make /s(0) smaller while still fitting lattice results at large q 2 . In figure |5| 
we show /s(0) as a function of the couplings g and g res for two choices of the extra resonance 
mass, m res = 6100 MeV, 8100 MeV. The problem in this case is that it is not possible to make 
a statistically acceptable fit to f£ and /g simultaneously. One could try to add an extra 
resonance in the J = channel also, but while our choice of m res = 6100 MeV for J = 1 may 



be motivated by potential models |H] or lattice results p2| , we do not know whether or how 
to set the mass for additional J = resonances, having already set the C values to account 
for rather poorly known resonances. This emphasises the importance of looking at / + and 
f° together, even though / + is the experimentally accessible form factor. 



5 Conclusion 

Our model is extremely simple, using only tree level HMChPT information for the two 
particle irreducible amplitude Vjj, thereby incorporating only the first excited hadron state. 
Furthermore we fix to a constant an allowed polynomial in q 2 muliplying the Omnes exponen- 
tial factor in equation (|J). Thus, deviations from LCSR results for / + are not unexpected 
because those calculations incorporate effects of higher resonances and continuum states. 
Taking our model beyond leading order is not possible at present because of the prolifera- 
tion of undetermined parameters which would appear in the next order of HMChPT and the 
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lack of experimental data to fix them. This is a standard difficulty in using effective theories 
at higher orders. 

The simple model presented here gives an excellent description of semileptonic .D-decays. 
For £?-decays it gives a good description of the lattice data near g^ ax and is also compatible 
within two standard deviations with LCSR predictions at q 2 = 0. Moreover, it provides 
a framework compatible with heavy quark symmetry, naturally accomodating pole-like be- 
haviour for / + and, simultaneously, non-constant behaviour for f°. Previously, as pointed 
out in 0, @], a difficulty for form factor models with pole- type behaviour for / + was fixing 
a behaviour for f° which satisfied both the relation /°(0) = / + (0) and the requirements 
of heavy quark symmetry. Pole-like behaviour of / + turns out again to be feasible in our 
model, thanks to the fact that the B* is a bound state rather than a irB resonance. 

Qualitatively, the results found here are encouraging. However, the larger value found for 
/s(0) compared to that from LCSR calculations would lead to appreciably smaller values 
for \V u b\- We caution the reader that this should not be taken to indicate a large theoretical 
spread in the value of \V u b\ from exclusive semileptonic B — > 7r decays: one should bear in 
mind the simplicity of the model used. We indicated how a second resonance in the J = 1 
channel can restore compatibility with both LCSR and lattice results for f£, although this 
shifts the problem to making fg compatible with the lattice data in a combined fit and 
emphasises the importance of using information from both form factors. 
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